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Given A'^ quantum systems prepared according to the same 
density operator p, we propose a measurement on the TV-fold 
system which approximately yields the spectrum of p. The 
projections of the proposed observable decompose the Hilbert 
space according to the irreducible representations of the per- 
mutations on points, and are labeled by Young frames, 
whose relative row lengths estimate the eigenvalues of p in 
decreasing order. We show convergence of these estimates in 
the limit oo, and that the probability for errors decreases 

exponentially with a rate we compute explicitly. 



I. INTRODUCTION 

The density operator of a quantum system describes 
the preparation of the system in all details relevant to 
statistical experiments. Like a classical probability distri- 
bution it cannot be measured on a single system, but can 
only be estimated on an ensemble sequence of identically 
prepared systems. In fact, if we could determine the den- 
sity operator on a single quantum system, we could com- 
bine the measurement with a device re-preparing several 
systems with the measured density operator, in contra- 
diction to the well-known No-Cloning Theorem . This 
points to a close connection between the problem of esti- 
mating the density operator and approximate cloning. In 
the case of inputs promised to be in a pure state the op- 
timal solutions to both problems are known , and it 
turns out that in a sense the limit of the cloning problem 
for output number M — > oo is equivalent to the estima- 
tion problem. The "optimal" cloning transformation was 
shown in this case to be quite insensitive to the figure of 
merit defining optimality ||^. 

In the case of mixed input states much less is known 
about the cloning problem. It is likely that in this case 
there may be different natural figures of merit leading 
to inequivalent "optimal" solutions. Even the classical 
version the problem is not trivial, and is related to the 
so-called bootstrap technique in classical statistics. 

The estimation problem certainly has many solutions. 
In fact, any procedure of determining the density ma- 
trix through the measurement of the expectations of a 
suitable "quorum" of observables 1^, such as in quan- 
tum state tomography Q is a solution. Other methods 
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include adaptive schemes [|j where the result of one mea- 
surement is used to select the next one. In all these cases, 
the estimate amounts to the measurement of an observ- 
able on the full input state p®^, which factorizes into 
one-site observables. What we are concerned with here 
is, as in the work of Vidal et. al. the search for im- 
proved estimates, admitting arbitrary observables on the 
A^-fold input system, including "entangled" ones. In con- 
trast to 1^ however we are not interested in estimators 
which are optimal for a more or less general figure of 
merit, but in the asymtotic behaviour if the number N 
of input systems goes to infinity (in this context see also 
the work of Gill and Massar |]lO|). 

When H = is the Hilbert space of a single system, 
the overall input density operator of the estimation prob- 
lem is /o®^, which lives on the N^^ tensor power Ti!^^ . 
This space has a natural orthogonal decomposition ac- 
cording to the irreducible representations of the permuta- 
tion group of N points, acting as the permutations of the 
tensor factors. Equivalently, this is the decomposition ac- 
cording to the irreducible representations of the unitary 
group on Ti (see below). It is well-known that this or- 
thogonal decomposition is labeled by Young frames, i.e., 
by the arrangements of N boxes into d rows of lengths 
Yi > 1^2 >■■■,> id > with = N. There is a 

striking similarity here with the spectra we want to esti- 
mate, which are given by sequences of the eigenvalues of 
p, say, ri > r2, ■ • • , > > 0, with = 1. The basic 

idea of this paper is to show that this is not a superficial 
similarity: measuring the Young frame (by an observable 
whose eigenprojections are the projections in the orthog- 
onal decomposition) is, in fact, a good estimate of the 
spectrum. More precisely, we show that the probability 
for the error \Ya/N — ra\ to be larger than a fixed e for 
some a decreases exponentially as TV — > cxo. 

Our basic technique is the Theory of Large Deviations 
pl| , in particular a result by Duffield on the Large 
Deviation properties of tensor powers of group represen- 
tations. This will allow us to compute the rate of expo- 
nential convergence explicitly. 



II. STATEMENT OF THE RESULT 

In order to state our result, explicitly, we need to re- 
call the decomposition theory for iV-fold tensor prod- 
ucts. Throughout, the one-particle space H will be the 
d-dimensional Hilbert space C*, with d < oo. Two group 
representations play a crucial role: firstly, the represen- 
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tation X ^ of the general linear group GL((i, C) 

and, secondly, the representation p i— > S'p of the permu- 
tations p S Sat on iV points, represented by permuting 
the tensor factors: 

S'p Vl (X" • • • V'Af = V'p-il ® • • •'0p-iAf • (1) 

The basic result is that these two representations are 
"commutants" of each other, i.e., any operator on 7^"^^ 
commuting with all X®^ is a linear combination of the 
S'p, and conversely. This leads to the decomposition 

U^'^^^TZy^Sy (2) 

Y 

X«^ = 0^y(X)®l (3) 

Y 

Sp = ^l®iTY{p) . (4) 

Y 

where Try : GL(d,C) ^ S(7^Y) and ^fy : ^ B{Sy) 
are irreducible representations, and the restriction of Try 
to unitary operators is unitary. The summation index 
Y runs over all Young frames with d rows and N boxes, 
as described in the introduction. We denote by Py the 
projection onto the corresponding summand in the above 
decomposition. 

Let us consider now the estimation problem. As al- 
ready discussed in the introduction, we are searching for 
an observable Ej^i describing a measurement on TV d-level 
systems, whose readouts are possible spectra of d-level 
density operators. The set of possible spectra will be 
denoted by 

E = {,seR^|x>0, = (5) 

where x > denotes the ordering relation on R"* given by 

s > : Sj > Sj+i for all j = 1, . . . , d — 1. (6) 

Technically, En must be a positive operator valued mea- 
sure on this set, assigning to each measurable subset 
A C E a positive operator ^^(A) G B{H^'^), whose 
expectation in any given state is interpreted as the prob- 
ability for the measurement to yield a result s £ A. 

The citerion for a good estimator En is that, for any 
one-particle density operator p, the value measured on a 
state is likely to be close to the true spectrum r e E 
of p, i.e., that the probability 

KNiA):^tr{EN{A)p®^) (7) 

is small when A is the complement of a small ball around 
r. Of course, we will look at this problem for large N. 
So our task is to find a whole sequence of observables 
En, N = 1,2,..., making error probabilities hke (|7|) go 
to zero as TV ^ oo. 

The search for efhcient estimation strategies En can 
be simplified greatly by symmetry arguments. To see 



this, consider a permutation p G Sn- If we insert the 
transformed estimator SpEn{A)S* into Equation (0) we 
see immediately that Kn{A) remains unchanged. Re- 
placing En{A) by the average Wr^ Ep^s^ SpEn{A)S; 
shows that we may assume [En{A), Sp] = for all per- 
mutations p, without loss of estimation quality. A sim- 
ilar argument together with the fact that the quality of 
the estimate is judged by some criterion not depending 
on the choice of a basis in Ti. shows that we may as- 
sume in addition that En (A) commutes with all unitaries 
U®'^. But this implies according to Eq. (||) and (|) that 
En (A) must be a function of the projection operators 
Py : n®^ 7^y (g) 5y defined at the beginning of this 
section. If we require in addition that each En {A) be a 
projection, which is suggestive for ruling out unnecessary 
fuzziness. En must be of the form 

En{A)^ J2 Py, (8) 

where sn is an arbitrary mapping assigning to each 
Young frame Y (with d rows and N boxes) an estimate 
sn{Y) G E. In other words, the estimation proceeds by 
first measuring the Young frame projections Py and then 
computing an estimate sn{Y) on the basis of the result 
Y. 
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FIG. 1. Probability distribution tr(p®'^Py) for d = 3, 
N = 120 and r = (0.6, 0.3, 0.1). The set E is the triangle with 
corners A = (1, 0, 0), B = (1/2, 1/2, 0), C = (1/3, 1/3, 1/3). 

The simplest choice is clearly to take the normalized 
Young frames themselves as the estimate, i.e., 

SNiY) = Y/N. (9) 

It turns out somewhat surprisingly that with this choice 
the En (A) from Eq. (^) form an asymptotically exact 
estimator. By this we mean that, for every p, the proba- 
bility measures Kn from Eq. (0) converge weakly to the 
point measure at the spectrum r of p. Explicitly, for each 
continuous function / on E we have 

hm J^f{s)KNids) = 

= ^^^T.f{^)Hp''''Py) = m- (10) 

We illustrate this in Figure ^ for d = 3, = 120, and 
p a density operator with spectrum r = (0.6,0.3,0.1). 
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Then S is a triangle with corners A = (1,0,0), B = 
(1/2,1/2,0), and C = (1/3,1/3,1/3), and we plot the 
probabilities tr(p®^Py) over Y/N e E. The explicit 
computation uses the Weyl character formula |]l3|, IX. 9.1] , 
which we do not need elsewhere in the paper. 

Clearly, the distribution is peaked at the true spectrum 
and our claim is that this will become exact in the limit 
N oo. To prove convergence we will use large devia- 
tion methods which give us not only the convergence just 
stated but an exponential error estimate of the form 



Kn( A) « exp 



-N inf I(s) 

sGA 



(11) 



where / denotes a positive function on E, called the rate 
function, which vanishes only for s = r. 

For the statement of the main theorem we say that a 
measurable subset A C S has "small boundary" , if its 
interior is dense in its closure. A typical choice for A is 
the complement of a ball around the true spectrum. 

Theorem. The estimator deSned in Eqs. and is 
asymptotically exact. Moreover, we have the error esti- 
mate 



lim —In i^Ar(A) = inf I(s), 



(12) 



for any set A C S with small boundary, where the rate 
function / : S ^ [0, oo] is 



/(.) 



(In Sj — In rj 



(13) 



The expression for / is the relative entropy ||lj| of the 
probability vectors s and r. Relative entropies occur also 
as the rate functions in Large Deviation properties of in- 
dependent identically distributed (classical ||T^ or quan- 
tum ||l^) random variables, although there seems to be 
no direct way to reduce the above Theorem to these stan- 
dard setups. 



III. SKETCH OF PROOF 

Rather than giving a proof of every detail, our aim here 
is to explain why the the scaled Young frames Y/N ap- 
pear in the estimation problem. The crucial observation 
is that the Young frame (Yi, . . . , Yd) is the highest weight 
of the representation Try in the ordering t> and this order- 
ing is directly related to picking out the fastest growing 
exponential in certain integrals of the measures K^. 

The integrals we need to study are the Laplace trans- 
forms of the measures K^. We introduce the "scaled 
cumulant generating function" 



c(ii) — lim — 



In 



KNids] 



(14) 



where rj G IV^, and 77 • s is the scalar product. If the 
measures Kn behave like Eq. (^l|) the integrand near s 
behaves like exp N{r] ■ s — I{s)), and the largest contribu- 
tion comes from the fastest growing exponential: 



c{rj) = sup(7/ • s — /(s)) 



(15) 



This is an instance of Varadhan's Theorem u%, which 
has a converse, the Gartner-Ellis Theorem ||1 1| , Thm. 
II.6.1]: if the limit ( p^ ) exists, and is differentiable then 
the estimate in the Theorem holds, with the rate func- 
tion determined from (|l^) by inverse Legendre transfor- 
mation. We will follow DufReld ||lj] by computing the 
limit (^4|) from group theoretical data. 

Consider the "maximally abelian subgroup" C C 
GL{d, C) of diagonal matrices 



Ph = diag(exp(/ii) 



,exp(/id)) 



(16) 



for h G C*. Since these commute, all the opera- 
tors TTyiPh) commute in every representation Try, and 
can hence be simultaneously diagonalized. The vectors 
fi = {fjLi, . . . , jjLd) such that TrY{ph)'4' = exp(/x • h)4' for 
some non-zero vector ip are called weights of the repre- 
sentation Try. The dimension m{p) of the corresponding 
eigenspace is called the multiplicity of p. One particular 
weight (with multiplicity one) is the Young frame Y itself 
(interpreted as an element of R'') and it turns out that 
Y is the maximum (the "heighest weight") among all 
weights of Try, in the >-ordering from Equation (p). Rep- 
resentation theory of semisimple Lie algebras [13 1 shows 
that each irreducible, analytic representation of GL{d, C) 
is uniquely characterized (up to unitary equivalence) by 
its highest weight Y. 

In order to estimate the integral (|lj), we need the 
quantities tr(p'^^Py). For simplicity we assume that p is 
non-singular, i.e., an element of GL(c?, C). By Equation 
(^ we have 

tr(p®^Py) = tr(7ry(p) l)= xy{p) dim(5y), (17) 
where 



Xy{p) ■■= tr(Try(p)) 



(18) 



is the character of the representation Try. Since xy is 
unitarily invariant {xy {U pU*) = Xy{p)) we may assume 
without loss of generality that p is diagonal and its matrix 
elements are arranged in descending order. Using the 
notation from Equation (|l^) this assumption reads 

p = Ph <^ C with /i > and ^ exp(/ij) = 1. (19) 



Hence we can express xy{p) in terms of the weights of 
Try: 



XY 



{p) ^^m{p)exp{p ■ h), 



(20) 
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where the sum is taken over ah weights of Try. Since 
hi> and Y t> fj, for ah ii we see that exp(F • h) is the 
largest exponential We therefore estimate 

exp(r • h) < xy{p) < dim (7^y) exp(y • h). (21) 

We will combine this with the consequence of Weyl's di- 
mension formula that dim(7?.y) is bounded above by a 
polynomial p{N) in TV, uniformly in Y Lemma 2.2]. 
Hence, for any h,ri E R"*, ft,, 77 > the two expressions 

= Y.XY{ph) e'T'' 6XTn{SY) (22) 

Y 

and 

J'(/i, r;) = ^ e^h+nVY dini(5y) (23) 

Y 

are asymptotically equivalent in the sense that 
(l/A^) (In J(/i, 77) — In J'{h,ri)) 0. In the same sense 
we can continue the chain of equivalences 

J{h, Tj) « J'{h, 77) = f{h + r;, 0) « J{h + ri,Q) = 

= Mp^^,)-(trp.+,)^. (24) 

Hence, if = cxjp{ha) are the eigenvalues of a non- 
singular density operator, we get for Eq. (|l^ the expres- 
sion 

c{r]) = In ^ Ta ex.p{ria) (25) 

a 

It is then a simple calculus exercise to verify the above 
rate function as the Legendre transform /(s) = sup^(ry ■ 
s - c{r])). 

This concludes our sketch of proof. In order to expand 
it into a full proof, one needs to extend the computation 
of 0(77) to rjf 0, and prove that this extension has the 
required regularity properties for the application of the 
Gartner-Ellis Theorem cited above. This has been car- 
ried out by Dufheld ||l^ in a context which is on the one 
hand wider, because it includes tensor powers of much 
more general representations of semisimple Lie groups, 
but on the other hand is narrower, because it contains 
only the case p = d^^l of our Theorem. However, one 
can extend Duffield's result by multiplying his measures 
Km with the factor Xv(p)/xy(l), and using for this fac- 
tor the estimate (^l|). 

IV. DISCUSSION 

Although the estimate we discuss is asymptotically ex- 
act, it is not at all clear whether and in what sense it 
might be optimal, even for finite N. We have experi- 
mented with various figures of merit for estimation and 



found different "optimal" estimators for low N, rarely 
coinciding with the En determined by (^). It is also not 
at all clear how much could be gained by optimization 
here. 

An interesting extension will also be the construction 
of estimators for the full density operator. It is very 
suggestive to compose this out of the above estimator 
for the spectrum, and to use for each Young frame a 
covariant observable to estimate the eigenbasis of p. The 
density of the covariant observable might be based on the 
highest weight vector of tty ■ 
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